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1. INTRODUCTION 

Generally, electric power system is studied in terms of generation, transmission and distribution systems 
in which all generators are operated synchronously at nominal frequency to meet the demand load. The frequency 
deviation in the power system is mainly due to mismatch between the generation and load plus losses at every 
second. There may be small or large frequency deviation based on the mismatch between generation and load 
demand. These mismatches due to random load fluctuations and due to large generator or power plant tripping out, 
faults etc. respectively. However, deviations could be positive or negative. The role of load frequency control(LFC) 
is to mitigate frequency perturbation. Thus the power system will operate normally [1, 2]. This can be achieved 
by adopting a auxiliary controller in addition to the primary control(Governor). From literature[1, 2], conventional 
controller is used as auxiliary or secondary control. To get the parameters of this controller, the power system is 
modeled and simulated using MATLAB. This paper deals with the modeling of power system through fractional 
order differential equations and design of controller. 

The fractional order dynamic system is characterized by differential equations in which the derivatives 
powers are any real or complex numbers. The approach of fractional order study is mainly used in the area of 
mathematics, control and physics [3]. The precision of modeling is accomplished using the theory of fractional 
calculus[4]. In view of above fact, integer operators of traditional control methods have been replaced by concept of 
fractional calculus[5, 6]. Many modern controllers for LFC as secondary controllers are available like sliding mode 
control[7], two degree of freedom PID controller[8], fuzzy controller[9], microprocessor based adaptive control 
strategy[10] and direct-indirect adaptive fuzzy controller technique[11, 12, 13]. It can be observed that power system 
parameters may alter due to aging, replacement of system units and modeling errors, as a consequent problem 
to design a optimum secondary controller becomes a challenging work. From literature, it is noticed that robust 
controller is inert to system parameter alteration. Thus, a good robust controller design is needed to take care of 
parameter uncertainties as well as load disturbance in power system. 

In literature, lot of robust control methods are presented for disturbance rejection and parameter alteration 
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for LFC. Still a vast research is going on internal model control(IMC) by researchers due to its simplicity. With the 
two degree of freedom IMC(TDF IMC)[14], both set point tracking and load disturbance rejection can be achieved. 
As a consequent, IMC controller design is an ideal choice for secondary controller for LFC. Saxena[15] designed a 
TDF IMC for LFC using approximation techniques like Pade’s and Routh’s, which motivated to adopt the fractional 
IMC-PID controller as a secondary controller and is design based on a fractional reduced order model of a system. 


2. REDUCTION METHOD FOR SINGLE AREA THERMAL POWER SYSTEM 
This section deals with framing of fractional order model of a single area power system using a step error 
minimization technique through genetic algorithm. This is segregated into following subsections. 


2.1. System investigated 


The proposed work deals with modeling of the power system to design secondary controller. Due to this pur- 
pose, a single area non-reheated thermal power system has been considered[1]. The thermal power system equipped 
with different units like generator G p(s), turbine G;(s), governor G,(s), boiler etc. and their dynamics are given 
by (1) 


(1) 


Figure 1. Single area power system linear model 


The block diagram of a single area power system is shown in Fig. 1, where AP, is Load disturbance(in 
p.u.MW), AXc is Change in governor valve position, A Pg is Change in generator output(in p.u.MW), u is Control 
input, R is Speed regulation(in Hz/p.u.MW) and A f(s) is Frequency deviation(in Hz). 

The overall transfer function is attained as 

Casel: From Fig. 1 assume Af(s) = Afi (s) when APi(s) = 0, the corresponding transfer function is G(s). 
Case2: From Fig. 1 assume A f(s) = A fo(s) when u(s) = 0, the corresponding transfer function is G2(s). 
Applying the theory of superposition principle to power system model, the overall transfer function is given by (2) 


Af(s) = Afi(s) + Afe(s) = Gi(s)u(s) + Go(s)APi(s) (2) 
The aim is to find control law u(s) = —K(s)Af(s) which mitigates the effect of load alteration on frequency 
deviation, where K (s) is fractional IMC-PID controller. 


2.2. Fractional system representation 


This subsection deals with the fractional order systems through which it can develop a proposed model for 
power plant to design a secondary controller. 
The representation for a linear time invariant fractional order dynamic system[16] is given as , 


H (D1 )y(t) a F(D%0P1---Bm u(t) (3) 
H(D%---om) — 5 aD% , F(DP081.---Bm) = 5 by DPF (4) 
ars k=0 
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where y(t) and u(t) are output and input vectors respectively and D is differential operator. a, Bp are the order of 
derivatives. ax and by, are coefficients of derivatives, ag, bgeR. Here H, F Fractional dynamic systems. 

The transfer function of fractional order dynamic system is obtained by applying Laplace transform to (3) and (4) 
(initial conditions are zero) and is given as (5) 


m ajk 
G(s) = k=0 by (s ) 


= 2ek=0 rls) 5 
Ene a 


2.3. Design of proposed system 


The full order transfer function of single area power system is obtained from (1) and (2), which is given by 


(6) 


Kp 
— TpTrTas? + (TpTr + TrTg + TaTp)s? + (Tp + Tr + Ta)s + (1+ Kp/R) 


Gi(s) (6) 


substitute the values of Tp = 20 sec, Tr = 0.3 sec, Tg = 0.08 sec, R = 2.4, Kp = 120 from [14] in (6), we get G4 (s) 
as (7) 


250 
s? + 15.885? + 42.465 + 106.2 


A(s) = Gi(s) = (7) 


The equation (7) represents integer higher order model which is converted to fractional order model assumed to be 
A'(s) 
Consider fractional SOPDT reduced model A’ (s) given by (8) 


Kye" 


A' (s) = ————_ 
(a sot psc +q 


(8) 
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Figure 2. Comparison of step responses of full order model with fractional SOPDT, Routh and Pade approximation 


Here the order of A’(s) is less than the A(s) and is in fractional form. The step error minimization 
technique[17, 18] through Genetic Algorithm(GA)[19, 20] is utilized to obtain parameters of A’(s) are given as 
Kı = 16.385, L = 0.095, b = 1.897, c = 0.962, p = 1.954, q = 7.031. This is done through MATLAB & simula- 
tion. 

Thus the attained reduced fractional model A’(s) is 


16:385e 0-09? 


A'(s) = 
(s) = -1897 4 195450962 1 7.031 


(9) 


The step responses of original model, proposed, Pade’s approximation and Routh’s approximation are compared and 
shown in Fig. 2. The performance index ITSE of proposed Pade’s and Routh,s methods are 0.0015, 0.027 and 0.06 
respectively. From Fig. 2 and above ITSE values, it is observed that the response of proposed method is very closer 
to full order model(original). 
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Figure 3. Block diagram of a) IMC configuration b) equivalent conventional feedback control 


3. FRACTIONAL IMC CONTROLLER DESIGN 
3.1. Internal Model Control 


In this section, model based IMC method for load frequency controller design is considered, which is 
developed by M. Morari and coworkers[21, 22]. The block diagram of IMC structure is shown in Fig. 3a, where 
Crmc(s) is the controller, P(s) is the power plant and P, (s) is the predictive plant. Fig. 3b shows block diagram 
of conventional closed loop control. From Fig. 3a and Fig. 3b, we can relate C(s) and Crarc(s) as 


Cruc(s) 


C(s) = (10) 
is) 1 — Crmols)Pm(s) 
Steps for IMC controller design[23] are as follows. 
Step1: The plant model can be represented as 
Pm = Pin Pin (11) 


where P7 = minimum phase system and P= non-minimum phase system, like zeros in right side of S-plane etc. 
Step2: The IMC controller is 


1 1 
— f(s), f(s) = (4 rast)" 


Cruc(s) = (12) 


where f(s) is low pass filter with steady state gain of one 

where 7, is the desired closed loop time constant and r is the positive integer, r >1, which are chosen such that 
Crmc(s) is physically realizable. Here r is taken as 1 for proper transfer function. 

The FIMC controller is designed for fractional SOPDT is given by (9) via method discussed below. 

Consider the system [23], is given by (13) 


ke 9s 
asf +a,s% +1 


where a : 0.5 < a < 1.5, 8: 1.5 < 6 < 2.5, 0= time delay. Here az = 7? and a, = 2r. 
Using (11), the invertible part of Pn (s) is 


Pm(s) = B>a (13) 


k 


Pals) = —r 14 
nl) = ae (14) 
Using (12), the Fractional IMC controller is 
B a 
azs” +a;s*+1 1 
G = 15 
tuc(s) ke ie) (15) 
substitute (15) in (10), then the conventional controller C(s) is evaluated and simplified as 
B a 
ags’ +ais% + 1 
C(s) = mm 16 
(s) k(t-st+1 + Os) (16) 
where e~°* is approximated as (1 — @s) using first order Taylor expansion[23]. 
Again C (s) is decomposed into FIMC PID filter via technique discussed below. 
Multiplying and dividing RHS of (16) by s~° 
B a 1 —a 
oe (ags’ +a1s* +1) 8 (17) 


k(Tesòtt + 8) s797 
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substitute az = T? and a; = 2Ċr in (17) and rearranged as, 


s172 2C7 1 T 
B-a 
T+ e/a he + sera + a¢8 


C(s)=[ )] (18) 


where first part is fractional filter and second part is fractional PID controller. 


3.2. FIMC to single area non-reheated power system 


To design FIMC, the (9) can be re-framed in the form of (13) is 


2.3303e 70-0958 


A'(s) = 19 
(9) = 974951897 40.2676502 4 1 a 
From (13),(18) and (19), we get 7 = 0.3768, 0 = 0.095, k=2.3303, 6 = 1.897, a = 0.962, ¢ = 0.3551 
substituting above values in (18), the fractional IMC-PID filter for single area power system is 
0-038 
C(s) 1.21(1 + 3.736s— °°? + 0.530559?) (20) 


~ 1+ 10.5267.5> 


The value of 7, and A are selected in such a way, that it minimizes tracking error and achieves robust performance. 


4. RESULTS AND DISCUSSIONS 

In this section, a model with proposed system and its controller is designed in simulation and an extensive 
simulation is carried out. In this model an performance index ISE is accompanied to determine the error of frequency 
deviation. Based on ISE, Overshoot, undershoot and settling time, we choose best À and 7,. The obtained parameters 
A and Te are A=0.22 and 7,=0.02 


Frequency Deviation(Hz) 


Figure 4. Comparison of response of power system using proposed with other methods 


To evaluate performance a step load disturbance A P,(s)=0.01 is applied to a single area power system 
as shown in Fig. 1. The frequency deviation A f(s) of proposed method in comparison with Pade’s and Routh’s 
method under nominal case is shown in Fig. 4. It is clear from Fig. 4 that the frequency deviation of the system for 
proposed controller due to load disturbance is diminished compared to Pade’s and Routh’s approximation methods. 
The performance index of proposed and other two methods are compared and shown in Table | under nominal case. 
From Table 1 it is observed that the performance index ISE is significantly low as compared with other methods. 


Table 1. Comparison of performance index for proposed and other reduced models under nominal and 50% Uncer- 
tainty cases 


Nominal case 50% Uncertainty case 
Methods Lower bound Upper bound 
ISE ISE ISE 
Pade’s approximation(Saxena) 8.4* 10-7 9.1* 10-3 8.4% 10° 
Routh’s approximation(Saxena) 8.7* 107 9.2 x 10% 8.9 x107 
Proposed method 1.410 ° | 8.44+10® 6.8% 10-° 
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4.1. Robustness 


Examining robustness of controller is vital because modeling of system dynamics is not perfect. So we have 
chosen fifty percentage uncertainty in system parameters to check robustness of controller. The uncertain parameter 
ði, for all ¿ = 1,2,..,5 are taken as[14]. Here ô is parameter uncertainty. The lower bounds and upper bounds of 
system uncertainty responses for proposed, Pade’s and Routh’s method are shown in Fig. 5a and Fig. 5b respectively. 


x10 x10 
a5 a5 
È 3 
E E 
5 a > Y. 
D E D We 
Q E Proposed A Il 27 Proposed 
zo zo W 

| E a a 's roximati C E a a 's approximati 

z -10 # Saxena Routh's approximation | | z -10 E Saxena Routh's approximation | | 
2 — —-—- Saxena Pade's approximation pi => >; Saxena Pade's approximation 
E 15 , 7 n n © 15 fl n 7 I 

0 2 4 6 8 10 0 2 4 6 8 10 

Time (s) Time (s) 


(a) (b) 


Figure 5. Comparison of response of power system using proposed method with other method for a) lower and b) 
upper bound uncertainties 


From Fig. 5a and Fig. 5b , it is noticed that the proposed controller is robust when there is plant mismatch 
and system parameter uncertainty. The performance index ISE of proposed and other methods under 50 % uncertainty 
case are compared and given in Table 1. It is observed that, there is significant difference in ISE between the proposed 
and Pade’s, Routh’s. Thus proposed controller is more robust compared to other two controllers. 


4.2. Proposed method extended to two area power system 


Areai 


+ y = 
u; I 1 + Ka Afi 
>| Cs >| > 7 ial F g 
ACE, if 1+ sTo; APG|1+ sT, AP,; 1+sT,, 
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D 
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Figure 6. Block diagram of multi area power system 


The block diagram of multi area power system linear model is shown in Fig. 6. The parameter values of 
model are given below[24]. 
Tp, =T p2 =20 secs, Tp, =Tr2= 0.3 secs, Tg, = Ta2= 0.08 secs, Ry = R= 2.4, K pı = Kp2=120 
Here i=1,2 and two area power system is assumed to be identical for simplicity. 
As followed above subsections, the controller is designed for two area power system with an assumption that there 
is no tie line exchange power(Ti2 = 0). 
The resulting FIMC-PID controller for areal and area2 of two area power system are given as (21) 


50-038 
~ T+0.0180-64 
To evaluate the performance of controller, a step load disturbance A P;(s)=0.01 is applied to a two area 
non reheated power system. The frequency deviations in areal & area2 and deviation in tie line of proposed method 


C1(s) = C2(s) 1.21(1 + 3.736s—° 9° + 0.5305s°9%°) (21) 
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Figure 7. Responses of two area power system 


are compared with Wen Tan method[24], which is shown in Fig. 7. It is clear from figures that the deviations of the 
power system for proposed controller due to load disturbance is diminished compared to Wen Tan method. 


5. CONCLUSION 

A good robust LFC technique is required to act against load perturbation, system parameter uncertainties 
and modeling error. In this paper a good approximation model reduction technique i.e step error minimization 
method is adopted to design a robust fractional IMC based PID controller for non-reheated thermal power system. It 
consists of fractional filter and fractional order PID. The tuning parameters, time constant 7, and non integer À are 
evaluated to get fast settling time and better overshoot/undershoot respectively. The simulation results showed that 
the proposed controller is more robust and good at set point tracking and for disturbance rejection. The performance 
of proposed method is good when applied to two area power system. 
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